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THE GENERALIZED RIEMANN PROBLEM FOR LINEAR 

DIFFERENTIAL EQUATIONS AND THE ALLIED 

PROBLEMS FOR LINEAR DIFFERENCE AND 

g-DIFFERENCE EQUATIONS. 

By George D. Birkhoff. 
Eeceived June 9, 1913. 

The program of obtaining a characterization of a function in sim- 
ple descriptive terms which are independent of the equations of 
definition of the function is a familiar one. To Riemann is due 
the formulation of this characterization for the algebraic functions 
and for the functions defined by ordinary linear differential equations 
without irregular singular points. In both of these instances the 
characterization involves a certain number of characteristic constants 
— the monodromic group constants in the last mentioned instance. 
Riemann also proposed the associated problem of assigning these 
constants at pleasure. 1 

During the last few years I have discovered that the program 
admits of extension in a number of directions. The aim of the present 
paper is to solve the generalized problem of Riemann for ordinary 
linear differential equations with irregular singular points, and the 
analogous problem for linear difference equations and for linear g— dif- 
ference equations. The formulation of the first and second of these 
problems has been given by me earlier. 2 At about the same time as 
myself, Norlund, in his fundamental work on linear difference equa- 
tions, was led to formulate essentially the second problem. 3 The 
third is stated in the present paper. 

The problem of Riemann for linear differential equations in its 

1 Werke, (zweite Auflage) pp. 37-39, 67-69. 

2 Trans. Am. Math. Soc, 10, 436-470 (1909), and 12, 243-284 (1911). These 
two papers will be referred to as I and II respectively. 

3 Memoires de l'Acad6mie Royale des Sciences et des Lettres de Danemark, 
series 7, 6, 309-326 (1911); C. R. vol. 156, pp. 200-202 (1913). 

In the second of these papers, Norlund gives a formulation and explicit 
solution of the hypergeometric difference equation problem. 
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classic form was first solved by Hilbert. 4 His treatment and Plemelj's 
elegant completion thereof 5 reposed alike upon a certain theorem 
whose proof was made by means of the Fredholm theory. Owing 
to the deep-seated analogy between linear differential and difference 
and g-difference equations, I have been able to apply a convenient 
extension of the same theorem in all cases; my proof is based on a 
method of successive approximations. 

Inasmuch as I have been able to simplify Hilbert's and Plemelj's 
treatment of the classic Riemann problem, I have ventured to include 
my treatment of it also. 

Part I. The Preliminary Theorem. 
§ 1. Some Definitions. 

Let C be a simply closed analytic curve in the complex x-plane. 
If the arc length along this curve from a fixed to a variable point is 
measured by s, and if I be the length of C, it is clear that a: is a single- 
valued analytic function of s with period / for s real, and that dx/ds 
is not zero. Consequently if we introduce a new variable r defined 

by 

2ir V— Is 

t = e l , 

a one-to-one analytic correspondence is set up between the points 
of the unit circle | r [ = 1 in the r-plane and the points of C. It will 
therefore be possible to choose p > 1 so that the circular ring in the 
r-plane, 

-<ItI < 



p ' 



= \ T = P» 



is transformed in a one-to-one and conformal manner into a ring in 
the rc-plane bounded by simply closed analytic curves C\ and C%, 
within and without C respectively, while at the same time the circle 
| r | = 1 is transformed into C. Let r = r (x) be the function which 
effects this transformation. 

Also let a (x) be any function continuous together with its deriva- 
tives of all orders along C, 6 and analytic save at a finite number of 

4 Gott. Naohr. (1905), pp. 307-338. 

5 Monatsh. f. Math. u. Phys., 19, 205-246 (1908). 

6 By definition we take df{x)/dx along a curve L as follows: 

<m = M/iHw (s>onL) . 

dx x ~ x x — x 
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points of C. These restric-, 
tions on a (x) ensure that we 
can choose regular curves D\ 
and D 2 , within and without 
C respectively, and osculating 
C at the points where a (x) is 
not analytic (Fig. 1) in such a 
way that on the continua lim- 
ited by D h D 2 , we have 




(1) 



a (x) \ £ K, 



a (x) — a (x) 



Fig. 1. 



^K; 



in these continua a (x) is defined as the analytic extension of a (x) on 
C. It is possible to extend further the definition of a (x) throughout 
the ring formed by C\, C 2 in such wise that inequalities of the type 
(1) hold; for this purpose it is clearly sufficient to choose real and 
imaginary components that join on continuously to the like com- 
ponents of a (x) along Di and D 2 , and to make each component satisfy 
inequalities of the same nature as (1). Such a choice can always be 
made. 

§ 2. On a First Type of Integral. 



Let us turn now to consider the integral 



(2) 



2ttV— l 



Jc t — x 



(t) 



a(t)dt, 



where g + (x) is a function analytic within C and continuous along C, 
and p is zero or a positive integer. Following Plemelj (loc. cit.) we 
shall term a function g + (x) of this description a regular inner function 
and affix to it a superscript + ; likewise a superscript - will indicate 
that a function is a regular outer function, i. e. is analytic in the 
extended plane without C, and continuous along C. 

We can demonstrate at once that the integral (2) represents a regu- 
lar inner function f+(x), or a regular outer function f-(x), according 
as x is within or without C. In the first place these functions are 
analytic within and without C respectively, as appears from (2). 

In the second place, by Cauchy's integral theorem we have 

(3) o--5== f^*^ L=f!M 4i 

2W— \Jc t—.x 2irV— lJci t — x 
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provided that x lies between C and C% since the function t p {x) g + {x) 
is analytic in the ring C, C\ and continuous along its boundary. Thus 
we may write 

(4) r {X) - 27V=TT X T (0 * (0 t=^ dt 

2tV— 1«/Ci t — x 

The first integrand on the right-hand side is continuous in x and t, for 
t on C and a; in the ring C, C 2 unless x = t, when the integrand is not 
defined; in the neighborhood of a pair of values x = t, the integrand 
remains finite by (1). Hence the first integral approaches a continu- 
ous limit as x approaches the boundary. Inasmuch as t is restricted 
to lie on d in the second integrand, the same statement is certainly 
true of the second integral. Hence f~(x) may be so defined as to be 
continuous along C. 

Likewise by means of the relation 



(5) T*MfM = ^hrJ c Tl 



dt 



* ■ i fifww* 

2tt V— \Jc x t — x 
valid for x between C and C\ by Cauchy's integral formula, we obtain 

(6) t (*) - r»(x)g+(x) a{x) = - J j== f ^ (t) g + (t) a( - t \~ a ^ dt 

2 7T "V 1 "C < * 



a{x) f t*M 
2ttVITiJci t 



t — X 



From this equation we can at once infer that f + {x) may be so defined 
as to be continuous along C. 

Thus / + (x) and f~(x) are respectively regular inner and outer func- 
tions. 

A comparison of the relations (4) and (6) which are both valid along 
C gives us the fundamental equation 

(7) / + (*) - /" 0) = t» (*) g + (*) a (x) along C. 

Let us now consider the maximum modulus of f~(x) outside of or 
along C. This maximum modulus, and likewise that for g + (x) within 
or along C, are attained on C of course. Suppose that we have along C 

(8) \g+{x)\^L. 
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Now modify the contour C of integration in (4) to Dj.. The integrand 
is analytic in / over the continua enclosed by C and D\, so that the 
value of the integral will not thereby be altered. (It must be remem- 
bered that x lies without C in (4).) From this modified form of (4) 
we obtain 



(9) l/-(*)|^™| f\7»(t)dt\+ flf^dt 



!■ 



upon applying (8) and (1). 

But the two integrals which appear in the right hand member of 
this inequality tend to zero as the unspecified integer p increases; in 
fact we have | r (a;) | < 1 within C so that t p (x) tends uniformly to 
zero in any closed continuum within C, as p becomes infinite. It is 
to be observed that the quantity [ t—x | which appears in the second 
integrand is never less than the minimum distance from C to C\, 
since x lies without C, and t is a point of Ci. 

These considerations demonstrate that for a given positive e, 
however small, the integer p may be chosen so large that for every 
regular inner function g + {x), we have 

(10) maximum of | f~ (x) \ ^ e { maximum of g + (x) } along C. 



§ 2. On a Second Analogous Type of Integral. 
In the same way we may treat an integral 

where g~{x) is a regular outer function, and p is zero or a positive 
integer. As before we denote the value of the integral for x within C 
by f + (x), and for x without C by f~(x). A discussion parallel to the 
earlier one in § 1 shows that f + (x) and f~(x) as thus defined are 
regular inner and outer functions respectively; in this case equation 
(4) is replaced by 



(4') f(x) = —J= (r-*{t)g-(t) 
2xV— \Jc 



a(t) —a(x) 

at 



t — x 
a(x 
27rV^TlJc2 t — x 
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and (6) likewise by 



(6') /" (aO+T-*(aO g-(x)a(x) = - ) f r^ (t) g' 

2t V— 1 Jc 



(t) a(t) t Z a x ix) dt 

, a(x)_ r r-" (t) r (*) dt 

2ttV— 1 Jc 2 t — x 

From these two equations there results at once 

(7') / + (3) - /" (x) = t-p (x) g- (x) a (x) along C. 

In order to develop an inequality for the modulus of / + (x) in this 
case, we note that the contour C in (4') may be modified to D%. The 
modification yields 

(90 \r(x)\^\ f\r-p(t)dt\+ rirlM^lj 

where L is the maximum of [ g~(x) j along C. But | t~ p (x) | tends to 
zero for x outside of C as p becomes infinite, since for such an x we 
have | t(x) \ > 1. We conclude therefore that for any positive e 
however small, the integer -p may be taken so large that for every 
regular outer function g~{x) we have 

(10') maximum of |/ + (x) | ^ e { maximum of \g~ (x) | } along C. 

A further property of the function /"(«), which is apparent from its 
definition, is that this function vanishes at x = oo . 

§ 3. Solution of a Pair of Matrix Equations. 

Throughout the present paper we shall be concerned with linear 
equations in n unknown functions, whose complete solution may be 
expressed in terms of n particular solutions. On this account we 
shall employ the matrix notation. 

We consider first a pair of matrix equations 

nu {F+(x)-F-(x) = T*>(x)G + (x)A(x), 

W \g-(x)- G+ (x) = t-v (x) F- (x) A-* (x) - I, al ° ng C " 

Here t(x) is the function defined in § 1; the matrix A(x) is a given 
matrix (a^x)) (i,j= 1, . . . ,n) of which each element is defined along 
C and has the properties specified in § 1 for the function a(x) (namely, 
it is continuous together with its derivatives of all orders along C, 
and analytic save at a finite number of points); furthermore the de- 
terminant | Aix) | is not to be zero along C. The symbol I stands 
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for the unit matrix (<5#) in which S a = 1, % = for i ?£ j, and A _1 {x) 
stands for the matrix inverse to A(x). The matrices F + (x), G + (x), 
F~(x), G'{x) are to be determined to satisfy (11), the first two as 
matrices of regular inner functions and the last two as matrices of 
regular outer functions. The matrix products on the right hand side 
are the customary matrix products, and the factors t p (x), t~ p (x) stand 
for the matrices (t p (x) 5y) and (r _p (x) 5 i3 -) respectively. 

It may be proved without difficulty that for p taken large enough 
a solution of these equations exists. To effect this proof we apply a 
method of successive approximations based on the sequence of 
equations 

Fr{x) = F ~(x) = G ~(x) = 0, G + (x) = I, 
i Fl+ (a .) _ Fl - ( x ) = t p ( x ) Go + ( x ) A ( x . ); 

(12) I Gr (x) - Gf (x) = t-" (x) F - (x) A- 1 (x) - /, * "' 

I Fr (x) - Fr (x) = r p (x) Or (x) A (x), 

Or (x) - Gr (x) = t-v (x) Fr (x) A' 1 (x) - 1, 



along C, 



The symbol is used to denote a matrix of zero elements, and the 
superscripts + and — are used to designate matrices of regular inner 
and outer functions respectively. 

If we write P -(x) = F ~(x) = 0, Q +(x) = G+(x) = I, and further- 
more 

(13) * Pm+ ( " T) = Fm+ (x) ~" Fm ~ 1+( - x ' ) ' Pm ~ ^ = F ™~ ^ ~ F <»-i~ ( - x ~)> 
} Q m + (x) = G m + (x) - G M _i + (x), Q m - (x) = G m - (x) - G m „r (x), 

it is clear that the sequence of equations (12) is equivalent to 



(14) 



P m + (x) - P m ~ (x) = r p (x) Q m _,+ (x) A (x), a i ong C , 

Q»T(x)-Q m +(x) = r-*(x) P m -r(x)A-\x), (m = 1, 2,'. . .). 

Here the superscripts are employed as before. 

The form of equations (14) is such that we can determine P m + (x), 
P m ~(x), Q m + {x), Q m ~(x) in terms of P m _r(x), Q m -i + (x) so that the 
mth pair of equations (14) is satisfied. In fact the first one of the 
rath pair of matrix equations may be broken up into ri l ordinary 
equations 

n 

PU,m + (x) — pij, n -(x) = t p (x) £ qi,\.m-i + (x)axj (x) along C, 

(i, j=l,.. .,«), 
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where the third subscript on the functions corresponds to the subscript 
on the matrix. We have already obtained a solution of an equation 
of the form 

/ + (x) — f~ (x) = t p (x) qi,\,m-i + (x) a X ; (a;) along C 

(compare with (7)) in the form of a definite integral. By forming 
the sum of f + (x) and f~(x) for X = 1, . . ,n we obtain for every i and j, 
elements Pij, m + (a;) and Pij, m ~ (x) which form the elements of P m + (x) 
and Pm~x( x ) w ith the desired property (14). 

Likewise we can break up the second one of the mth pair of matrix 
equations (14) into n 2 equations. A solution may here be built up in a 
similar way (compare (7')- It must be observed that since the de- 
terminant of A (x) does not vanish along C, the elements of A~ x {x) 
satisfy the conditions imposed on a(x) at the outset. 

Now if we recall the method of solution of (14), it is clear from 
(10), (10') that along C the maximum modulus of any element of 
P m ~ (x) or Q m + (x) does not exceed weZi m _i where L m _i denotes the 
maximum modulus of any element of P m _r(x) or Q m -i + (x) along C, 
and e is arbitrarily small uniformly for all values of m. This relation 
may be expressed in the simpler form 

(15) Lm ^ nelrm-i. 

The series formed by the elements in any ith row and jth column of the 
series of matrices 

P.- (ar) + Pi" (*) + .. ., Qo + (a-) + Qi + (x) + . . ., 

will therefore converge absolutely and uniformly provided that e is 
taken so small that ne < 1. 

But the sums of to + 1 terms of these two series of matrices are 
F m ~(x) and G m + (x) respectively, whose elements therefore converge 
uniformly to the elements of matrices F~(x) and G + (x) of regulafouter 
and inner functions respectively. If we recall that P ~(x) = 0, and 
that the integral form of representation of each element of P m ~ (x) 
(see (2)) makes each element of this matrix reduce to zero at x =' oo , 
it is plain that at infinity F~(x) reduces to the matrix 0. 

If we turn now to consider F m + ix) and G m ~ (x) along C, we see from 
what precedes and from equations (12), that these matrices also 
converge uniformly along C, and therefore respectively within C and 
without C, to matrices F + (x) and G~(x) of regular outer and inner 
functions. Since G ~(x) = 0, it is clear that G~{x), as well as F~{x), 
reduces to at x = °o . 
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The matrices F~(x), F + (x), G~(x), G + (x) thus obtained will satisfy 
{11), as appears from (12) by letting n become infinite. 

§ 4. Application to the Solution of a Single Matrix Equation. 

Multiply the second matrix equation (11) on the right by 
t v (x)A{x) and subtract it, member for member, from the first equation 
(11). There results 

(16) F+ (x) = t*> (x) [I + G~ (x)] A (x) along C. 

Inasmuch as G~(x) reduces to a matrix of zero elements at x = oo , 
the determinant of / + G~(x) and also of F + (x) is not identically zero. 
The matrix equation (16) admits of further simplification. In fact 
the function log r(x) is analytic along C and increases by 2xV — 1 
as x makes a positive circuit of C. If c lies within C the function 

<p (x) = log t (x) — log (x-c) 

is accordingly single-valued as well as analytic along C. But <f>(x) is 
of the form 

r" (x) g+ (x) a (x) (p = 0, g (x) = 1, a (x) = <f> (x)), 
so that by (7) we can find 6 + (x) and d~ (x) such that 

0+ (.t) — 6- (x) = (x) along C; 
this gives us 

t p (x) = (x-c) p e v <> + We pe -(x). 
Now let us write 

$ (x) = e -v»+{x) F+ (x), *(x) = (x-c) p e pe '^ [I + G~ (x)]. 

By these equations we define $(x) as a matrix of regular inner func- 
tions, and ^r(x) as a matrix of functions analytic outside of C except 
for a pole of order p at x = oo , and continuous along C; the determi- 
nant of neither $(x) nor \l>(x) vanishes identically. Between $(x) and 
^(x), by (16), we have the matrix relation 

(17) $ (x) = * (x) A (x) along C. 

It is this type of matrix equation which is important for the present 
paper. 
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§ 5. Further Properties of &(x) and ty(x). 

It is necessary for us to investigate further the nature of any solu- 
tion $ (x), ^ (x) of (17) in the neighborhood of the curve C. In the 
first place, it is to be observed that at points of G where all of the 
elements of A (x) are analytic, the same is true of the elements of "l'(x) 
and ^(x) ; in truth, the equation (17) shows us that analytic extension 
is possible across the curve at such points. 

We shall prove that the elements of these matrices possess derivatives 
of all orders, continuous at all points of C. 

Since the elements of A(x) have line derivatives of all orders along 
C we may write, for t and y upon C, 

(18) A(t) = A(y)+(t-y)~A(y)+... 

+ -^ |U (y) + (f- yfB (t, y), 

where the elements of B(t, y) are continuous functions of t and y 
along C. 7 Also by Cauchy's integral formula in matrix form we 
have from (17) for x within C 



cl k ~^ (x) = (— l) *- 1 
dx k - } 2 



W=T Jc(t-x)" M 



(k = 2, 3, . . .). 



If we substitute the above expression for A (t) in this last equation, we 
obtain a number of terms of the form (save for a constant multiplier) 

The integral is not altered in value if C is replaced by C 2 which lies 
outside of C. Therefore each of these terms represents a function 
analytic in x and continuous in y along C. There remains a single 
term not of the form (19), namely 



x V— 1 Jc\t — xj 



(- D *- 1 (ft - 

2 



7 When a differentiation or integration sign appears before a matrix, it is 
understood to apply to each separate element of the matrix. The stated 
property of B(t, y) comes at once from the explicit formula 



°«.»-M'(^iY d ^<- 
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If now x be made to approach a point Xq of C, and if y be taken as the 
foot of the normal from x to C, this term approaches a limit which is 
continuous along C. In fact the factor (t — y) k /{t — x) k remains finite 
for t along C, and approaches the limit 1 uniformly save in the vicinity 
of x$. Thus if x (and y also of course) approaches a-o, d k ~ 1 ^? {x)/dx k ~ l 
approaches a limit along C continuous for an arbitrary .r . 

A similar proof shows that ty(x) has derivatives of all orders, con- 
tinuous outside of and along C. This proof is based on the fact that 
the elements of A~ l (x) satisfy the restrictions placed on the function 
a{x) in § 1. 

The above results also lead to the conclusion that at any point 7 of C 
at which one or more elements of A(x) fails to be analytic, the ele- 
ments of $(x) and Sfr (x) admit of asymptotic expansion in a series in 
positive integral powers of x — y. This is an immediate consequence 
of an expansion like (18) for $(x) or ^(x) in which now t and y can 
be any points within or without C respectively, and B it, y) is contin- 
uous in £and y. 8 

§ 6. A Normal Form for $(«) and ty(x). 

By a series of simple normalizations it is always possible to obtain 
a solution $(.t), ^(x) of (17) such that$(x) is (as before) a matrix of 
regular inner functions, and ty(x) is a matrix of function analytic, 
without C in the extended plane except for a possible pole at x = 00 , 
and furthermore such that | $(3:) | does not vanish within or on C, and 
I ^(x) I does not vanish without C. 

This solution may be directly obtained from that found in § 4. 
If I $(x) I vanishes at x = c within C say, we can determine a matrix M 
of constants such that all the elements of the first row of M $(x) 
vanish at x = c, while \M\ 5^ 0. Now M$ (x), M^ (x) yield a new 
solution of (17), which has the properties given for $(.r), ty(x) in § § 4, 5. 
If we divide the elements of the first row of M <£(«) by x — c, we obtain 
a matrix $(a:) of functions analytic within C and continuous along C; 
if the same operation be applied to ^(a;), we obtain ^r(x), a matrix 
of functions analytic without C in the extended plane save for a pos- 
sible pole at x = 00 . Moreover $(.t) and ^(x) yield a solution of (17), 
since the effect of this operation is to alter the matrix on either side 
only in the removal of a factor x — c from the first row. By this 
device we have diminished the multiplicity of the zero of [ $(x) [ at 

8 For the relation between the existence of derivatives and of asymptotic 
series see W. B. Ford, Bull. Soc. Math. France, 39, 347-352 (1911). 
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x = c by one unit, without the introduction of further zeros of [ <£(*') | 
or | ty(x) | in the finite plane. 

An entirely similar process eliminates a zero of [ ^(x) | without C, 
or a zero of | $(x) | and | ^(x) | along C. In consequence of the 
results of § 5, if either of these functions vanishes along C, the other 
does also, both at least to the first order. 

This process may be continued so long as there remain zeros of 
<J?(x) or ^(x). It must however finally come to an end. If this were 
not the fact it would follow at once that both | <3?(x) | and | ^(x) | have 
a zero of infinite multiplicity at some one point of C where an element 
of A (x) fails to be analytic. But this cannot be the case, for let $(x) 
and M'(.r) be a solution of the following matrix equation 

(20) <§ (x) = A-'(x) ■§■ (x) along C, 

where the elements of 3>(x) and ^(x) are restricted like 3>(x) and ^(x) 
were found to be in § 4. The existence of such a solution becomes 
manifest by a mere interchange of the role of rows and columns in 
what precedes. Now from (17) and (20) we conclude that 

! $ (x) | • | $ (x) | = I * (x) I ■ I * (x) I along C. 

The function represented by either side of this equality is not identi- 
cally zero; and it appears from the two representations that it is 
analytic in the finite plane, and analytic or with a pole at infinity. 
Hence this function is a polynomial, and the multiplicity of the zeros 
of either | $(x) | or | ^f(x) | at any point of C is finite. 

When the process comes to an end the following result has been 
obtained: if the elements of A(x) are unlimitedly differentiate along 
C, analytic save at a finite number of points of C, and if | A (,x) | is 
not zero along C, there exists a solution $(x), ^(x) of the equation (17) 

$ (x) = ■* (x) A (x) along C, 

in which the elements of $(.r) are analytic within C, unlimitedly differ- 
entiable along C, and $(x) is of determinant not zero within or on C; 
and in which the elements of ^(x) are analytic without C in the ex- 
tended plane save for a possible pole at x = oo , unlimitedly differen- 
tiable along C, and ^(x) is of determinant not zero without C. 

Here the point x = co appears as an exceptional point. It is evi- 
dent that an arbitrary point x = a not on C may be used to take the 
role of the point at infinity. In fact a may also be taken to be a point 
of C. When this is the case, the elements of $(x) and ^(x) are finite, 
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or become infinite to finite order at x = a. To obtain such a solution 
&(x), ~$?(x) it is only necessary to divide each element of $(*), ty(x) 
by a suitable power of x — a, so chosen as to make each element of 
these matrices analytic at x = oo ; and then to apply the normaliza- 
tion above indicated, letting l/x — a replace x. The curve C can also 
be taken to be a simply closed analytic curve which passes through 
x= oo . 

The main part of the conclusion that has been deduced above was 
obtained by Hilbert and Plemelj (loc. cit.) with the aid of the Fredholm 
theory of linear integral equations. Independently of their work, I 
treated a special case (see I, § 1) which arose in a different form in 
connection with my study of the irregular singular points of ordinary 
linear differential equations. 

My proof in this special case suggested to me the above treatment 
of the general case by the method of successive approximations. 
The restrictions here placed on the elements of A (x) and on the curve 
C are not essential to this treatment, and I have very little doubt that 
these may be replaced by the weaker restrictions of Hilbert and Plemelj. 
Nevertheless I have been content to use the simplest restrictions 
consistent with the applications in view. 

A second proof in the special case has recently been given by me, 
Math. Ann., vol. 74 (1913) pp. 122 (see also Bull. Am. Math. Soc, 
vol. 18, 1911, p. 64). This second proof, which suggested itself to me 
at about the same time as the first, is practically the same as that given 
by Hilbert and Plemelj. To my considerable regret this relation- 
ship escaped my observation until it was too late for me to make 
suitable reference. 

§ 7. The Preliminary Theorem. 

The following is an extension of the preceding results which is 
convenient for the applications : 

Preliminary Theorem. Let C h . . ., C r be r simply closed analytic 
curves in the extended complex plane. Let Ai(x), . . ., A r (x) be matrices 
of functions defined and unlimitedly diffcrentiable along Ci, . . . , C r 
respectively, analytic save at a finite number of points of these curves and 
of determinant not zero. If furthermore at any point of intersection of 
C a , Cp, the matrices A a (x), A${x) are such that the formal derivatives of 
all orders of the matrix 

(21) A a (x)A fi (x)-A (x)A a (x) 
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vanish, there exists a matrix $(x) with the folloiving properties: 

(1) each element of <£(«) is analytic except along Ci,..., C r and at an 
arbitrary point x = a where the elements may become infinite to finite 
order; | $(a-) | nowhere vanishes save possibly at x = a; 

(2) the elements of $(x) are continuous and unlimitedly differentiable 
along each curve Ci from either side, analytic from either side save at 
points of intersection of the curves, or at those points where an element of 
Ai(x) fails to be analytic, or at x = a; if a lies on a curve C,-, the matrix 
(x — a) 1 Ai(x) [or x~ l A f (x) if a = <x>] is unlimitedly differentiable along 
C t for a suitable I. 9 

(3) if a + and — side of each curve C t are chosen, then 

JQ $ (x) = lJt? r $ (x)} A t (x t ) (i = 1, . . .,r), 

where the approach to the arbitrary point x t of C; is along the + and — 
side respectively. 

Let us begin by establishing the theorem in the case when a is 
not a point of C\, . . . , C r . It has already been established for r = 1 
(see § 6), with the single notational difference that $(.t) was replaced 
by either of two matrices $(a;) and \& (x), according as x was within or 
without C. To establish the theorem then, we need only show that 
if it is true for r ^ k it is also true for r = k + 1, when the theorem 
follows by induction. 

Assume that $^ (x) is the solution for r = k, and for the matrices 
&,..., C k , A-iix), . . . , A k (x), where C h . . ., Cj+i, and A x (x), . . ., A m (x) 
satisfy the requirements of the theorem for r = k + 1. Let us sup- 
pose for the moment that a solution $^ +1 (x) exists with the desired 
properties. If we write 

(22) **rt(aO = C(*)*»(aO, 

the following facts are clear from (1), (2), (3) of the theorem: (1') 
each element U (x) is analytic except along C\, . . . , C^i and at the 
specified point a, where its elements may become infinite to finite 
order; | C7" (a;) | nowhere vanishes save possibly at x = a; (2') the 
elements of U(x) are continuous and unlimitedly differentiable along 
each curve C ; from either side, analytic save at points of intersection 

9 A function will be termed unlimitedly differentiable at x = <x if when we 
write x = 1/x', the function of x' obtained by the substitution is unlimitedly 
differentiable at x' = 0. 
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of the curves or at those points where an element of A t (x) fails to be 
analytic; 

, ,. xix 4 + U(x) = j™- U(x) (i = 1, . . ., k), 

(y ) lim T T / n , / \ r lim 



,+ 



U(x)$k(x) = [ji!,,,- U(x)$k(x)}A k¥ i(x k¥ i). 



x=x k+i^ " v»> **v»' i x = x k+i 

The condition (3') necessitates that U (x) is analytic at any point 
of Ci (i = 1, . . . , k) and throughout the extended plane save at a and 
at points of C k +i- The condition (3') gives us in addition 

x=x k + 1 + U(x) = [,i,™- V (x)]$k(x k +i) Ak+i(x k+ i) $4"' (Xk+l). 

Conversely if U(x) satisfies these conditions (1'), (2'), (3') it is ap- 
parent that$ t+ i (a;), given by (22), will satisfy the conditions prescribed 
for $(«) in the theorem. 

But in view of the above simplification of (3'), the conditions (l')> 
(2'), (3') on U(x) are precisely those of the theorem on $(x) if we take 
r = 1, C = C k+ i, Ai{x) = $ k (x) A k¥ i{x) $k~ l (x)- To complete a 
proof we need only show that this matrix fulfils the conditions pre- 
scribed in the theorem along Ck+i; for then a matrix U(x) will exist 
which satisfies these conditions. 

The elements of the matrix 

(23) $k(x)Ak + i(x)$ k - l (x) 

are analytic at points of Cjt+i which are not points of intersection or 
contact with C h . . ., C k , or points at which an element of A i+ i(x) fails 
to be analytic; at these latter points the elements of $j(a;) are analytic 
so that the above matrix is unlimitedly differentiable in the neighbor- 
hood of such a point. The determinant of this matrix is equal to 
| Ak+i(x) | and nowhere vanishes along Ck+i. 

It remains only to examine the elements of the matrix near points 
of intersection or contact of C i+1 with one of the curves C\, ■ ■ ■ , C k . 
Suppose first that C k+ i intersects a single curve C { at such a point. As 
x moves along the curve C k +i and passes from the positive to the 
negative side of C,-, the matrix f> 4 (a;) changes to $t (x) A t (x), so that 
the matrix (23) changes to 

$* (a;) A ( (x) A k+1 (x) Ac 1 (ar) Qf 1 (x). 

Bearing in mind the condition of permutability imposed on 
Ai (x),. . ., A r (x) at such a point of intersection (see (21)), it becomes 
apparent that the elements of the matrix (23) are continuous at this 
point, and have equal backward and forward derivatives of all orders 
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at the point. Thus the elements of (23) are unlimitedly differentiable 
in the neighborhood of such a point along C i+1 . The same is true at 
more complicated points of intersection or points of contact, as a 
similar argument shows. 

The theorem is now demonstrated, at least for the case when a 
does not lie on a curve &,..., C r . If a does lie on such a curve, we 
first choose an a' which is not on these curves to replace a, and then 
resort to the simple device used at the close of § 6 to replace a' by a. 

Part II: The Problem of Riemann and its Generalization. 

§ 8. On Cauchy Matrices. 

Let T be a matrix of constants of determinant not zero, correspond- 
ing to the coefficients of a linear transformation. According to the 
well-known theory of classification of such matrices, based on the 
Cayley-Weierstrass elementary divisor theory, we may write 

T = C~ l I'C, 

where C is a matrix of constants, and in general I' is of the form (p ; %). 

Consider along with the matrix /', the matrix of functions 

I' (x) = 0A'5y) where 2xfc,V— 1 = logp, ( j = 1, . . ., n). When x 

makes a positive circuit of x = 0, I'(x) becomes I'(x)I'. The matrix 

T(x) = I'(x)C 

will accordingly alter to 

I'(x) I'C = I (x) C-C-U'C = T (x) T, 

when x makes such a circuit, i. e. T{x) will be affected by the prescribed 
linear transformation. 

Such a matrix Tix) or a simple modification thereof exists for every 
transformation T and is called a Cauchy matrix. 10 A certain measure 
of arbitrariness enters into the determination of T{x) when T is given; 
in general ki, . . . , k n are undetermined up to an additive integer. 
We note that the determinant of T{x) is not zero for x 9^ 0, 00 . 

A final property of the Cauchy matrices which is important for us 
is that they form the matrix solution of a linear differential system 

x^-LTix), 
ax 

10 Cf. Schlesinger, Vorlesungen uber lineare Differentialgleichungen, pp. 
122-140, where a complete treatment is given. 
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where L is a matrix of constants. In fact we have 
dT{x) _ dl'jx) n 

X ~ 7 — X 7 C 

ax ax 

But the matrix x dl'(x) / dx may be written as 

KI'(x), K= (kjSij). 
Thus we obtain 

x ^M = KI' (x) G = KG- C-W (x) C = LT (x). 

§ 9. The Monodromic Group Problem. 

The most elementary existence theorems for ordinary linear differ- 
ential equations show that the linear differential system 

(24) ~^- = R(x) Y(x) 

admits a. matrix solution Y(x) whose elements are analytic in the 
finite plane at every ordinary point where the elements of R(x) are 
analytic, and analytic at infinity if the elements of R(x) vanish to at 
least the second order at the ordinary point infinity; furthermore 
| Y(x) | is not zero at such points. All other points of the plane are 
called singular points in contradistinction to the ordinary points. A 
finite singular point at which the elements of R(x) are analytic or 
have a pole of the first order, or the point x = oo if each element of 
R(x) vanishes to the first but not always to the second order at 
that point, is termed a regular singular point. 

Suppose now that the elements of R(x) are rational, and that all of 
the singular points «i, . . . , a m are regular. These will be taken to lie 
in the finite plane. It is easy to show that the elements of Y(x) be- 
come infinite to only a finite order at x = eti, . . . ,a m . 11 When x makes a 
positive circuit of one of these points, Y{x) changes to Y (x) T it where 
T{ is a matrix of constants and | T t \ 9^ 0; in fact the most general 
solution is of the form Y(x) T where T is an arbitrary matrix of con- 
stants for which ] T \ 5^ 0, and after the circuit is made in the x-plane, 
Y(x) is still a solution of (24). 

If we start from a point x = c and make a circuit of a\, . .,a m in 



11 Cf. Trans. Am. Math. Soc, 11, 199-202 (1910). 
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such wise that the combined circuit is reducible to a point, it is clear 

that 

(25) T m T m -i ... Ti = I, 

a necessary relation between the matrices T\, . . . , T m . 

The problem of Riemann is the following: For assigned points 
ai, . . . , a m and assigned matrices T h ..., T m for which (25) holds, to 
construct a matrix Y(x) of functions of determinant not identically 
zero, with elements analytic save at <ii, . . .,a m where these elements 
are finite or become infinite to finite order, and undergoing a transfor- 
mation to Y(x) Ti as x makes a positive circuit of a { (i = 1, . . ., m). 

A solution of this problem has 
been given by Hilbert and com- 
pleted by Plemelj (loc. cit.). It is 
possible to obtain a more simple 
solution on the basis of the pre- 
liminary theorem. 

Let us surround «i, . . . , a m 
by small non-overlapping simply 
closed analytic curves Ci, . . . , C m 
and let us pass through Oi, . . . ,a m in 
cyclical order another closed ana- 
lytic curve D which meets each 
curve C t only twice (Fig. 2), in \ and m> say. 

Choose matrices A\. . .,A m , A m+ i= Ai of constants such that 




Fig. 2. 



(26) 



(i = 1, . . . ,m). 



Here A\, for example, may be taken at pleasure and A it . . . , A m are 
then determined. Define a matrix A{x) along D to be equal to Ai 
on that part of the curve which lies between C t and C i+1 [C m+1 = CJ 
and equal to 



(27) 



• rrn 



■ m 



Ai., + 



■k 



mi — U 



A; (i = 1, . . . ,m) 



for x on the part of Z> within C t . We will suppose that l u mi and D 
were so chosen that D does not pass through one of the finite number 
of point for which the determinant of any matrix (27) vanishes. 

The matrix A(x) as here defined is continuous along D, analytic 
save at the points where D intersects C h . . ., C m , and of determinant 
not zero along D. 
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By slightly modifying each element of A(x) along a small segment 
near either end of the arc of D within C t for i = 1, . . . , m, it is possible 
to obtain a matrix A(x) whose elements are analytic save at the end 
points of these segments, and unlimitedly differentiable along D. 
Such a matrix A{x) satisfies all of the conditions necessary for the 
application of the preliminary theorem in the case r = 1. 

Taking D as the curve and A(x) as the matrix, we can affirm the 
existence of a certain matrix $(x), possessing in particular the property 

(28) Jlf i+ $ (x) = [J%- $ (*)] A (a*), 

the approach to the point x\ of D being from within and without D 
respectively. Let us choose a — a m . 

If we extend $ (x) analytically across D between C,- and C i+ i, 
it becomes $(x)A t by (28). If we extend this matrix analytically 
back across D between C i+1 and C i+2 [0 m+2 = C 2 ], it becomes 
[$ (x) A i+ \~ l ] Ac 1 , or $(x) T t . That is, the matrix U (x) obtained 
from $(x) by analytic extension is analytic outside of Ci, . . .,C m and 
undergoes a transformation to U(x) T t when x makes a positive cir- 
cuit of a t . Furthermore the determinant of U(x) is not zero outside 
of Ci, . . . ,C m . 

Denote by Z^x—ai) the Cauchy matrix belonging to the trans- 
formation T t , so that Z i (x — di) undergoes a transformation to 
Z t (x — a,) Ti as x makes a positive circuit of a,-. Write 

(29) F(x) = Z(x)U(x), 

where F(x) is the solution of the Riemann problem to be constructed. 

The elements of Zix) must in the first place be single-valued and 
analytic outside of C\, . . . , C m , since U(x) undergoes the same trans- 
formation as that prescribed for F(x) about the points a h . . .,a m , and 
|C7(x)|^0. 

Furthermore within C,-, the elements of F(x) are to be analytic ex- 
cept at a ; where they may become infinite to finite order. Hence the 
elements of the matrix Y{x)Z i ~ 1 {x — a t ) must be single-valued and 
analytic within C,-, by the definition of Z { (x— a,), except for a possi- 
ble pole at x = a { . Along C t this matrix may be written 

(30) Zt (x) = Z (x) [U{x) Zi- 1 (x — m)] (i = 1, 2, . . . m). 

If we write $ (x) = Z(x) outside of C\,. . .,C m and also $(x) = Z,-(x) 
within C for i — 1,. . .,m, the equations (30) may be written 
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(31) x % $ (*) = [ x % $ (x)] A t (xd, A t (x) = U (x) Zci (x-ai). 

(i = 1, . . ., m). 

This suggests another application of the preliminary theorem, since 
the curves Ci, . . .,C m and the matrices A\(x), . . ., A m (x) of known 
functions satisfy the necessary restrictions. 

Let $(z) be the solution given by the theorem for a = a m , and let 
Zi(x),. . .,Z m (x), Z(x), be defined as equal to <£(«) within C\, . . ., C m 
and outside of these curves respectively. These functions will then 
satisfy (30). Let Y(x) be defined by (29). This matrix is clearly 
composed of elements analytic without and along Ci,..., C m , as are 
those of U(x); within C t , Y{x) continues analytically into 
Zi(x)Zi(x — Of) by (30) for i — 1,. . .,m, and consequently its ele- 
ments are analytic throughout the plane except possibly at ay . . .,a m ; 
its elements become infinite only to a finite order at a t since the 
elements of Z t {x — a t ) become infinite only to finite order at a ( . 
Furthermore by (29) Y(x) undergoes a linear transformation to 
Y (x) T { as x makes a positive circuit of a^. Thus the Riemann prob- 
lem has been solved. 

It is worthy of note that | Y (x) \ does not vanish for x 9^ o,- 
(i = 1, . . ., m). This is an immediate consequence of the fact that 
| Z\(x) \, . . ., | Z m (x) |, | Z(x) | do not vanish in their regions of defini- 
tions save at these points, and of the fact that the Cauchy matrix 
Z t (x — a,i) has a determinant which does not vanish save possibly at 
x = a t and x = oo . 

§ 10. A Generalization. Equivalence. 

A more general result can be deduced exactly as the results of § 9 
were obtained. Let us say that two matrices of functions Yi (x) and 
Yi{x) whose elements are analytic in the vicinity of * = a, but not 
in general single- valued or analytic at x = a, are properly equivalent 
at x = a if we have 

Y,{x) = A{x)Y,{x). 

where A(x) is composed of elements single-valued and analytic at x= a, 
of determinant not zero there; if this condition is not satisfied, 
but if the elements of A(x) have a pole or are analytic at x = a, let 
us say that Yi(x) and Y 2 (x) are improperly equivalent atx = a. 

This definition is convenient for the statement of the following 
result : Let ai, . . .,a m bem given points; let T\, . . ., T m be matrices of con- 
stants such that T m T m -i . . . Ti = I; let Zi (x), ..., Z m (x) be matrices of 
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functions analytic of determinant not zero in the vicinity of a t and under- 
going a transformation to Zi(x) T h . . ., Z m (x) T m as x makes a ■positive 
circuit of ai, . . ., a m respectively. There exists then a matrix Y(x) of 
functions of determinant not zero for x 5^ a\, . . . , a m and analytic 
save at these points, which undergoes a transformation to Y (x) T ( as x 
makes a positive circuit of a { (i = 1,. . .m); furthermore Y(x) is prop- 
erly equivalent to Zi (x), . . . ,Z m -i(x) at a h ... , o m _i and properly or im- 
properly equivalent to Z m (x) at a m . 

The result above stated is obtained when we replace the Cauchy 
matrices Z x ( x— ai), . . . , Z m (x— a m ) of § 9 by matrices Zi{x), . , . , Z m (x) 
having the properties specified. The line of attack is identical with 
that given in § 9. The facts concerning equivalence follow at once 
from the relations analogous to (30) : 

(32) Y (x) = Z { (x) Z t (x) (i = 1, . . . , m). 

Here Z t (x) is composed of elements analytic at x = a ; (i — 1, . . . ,m — 1), 
and of determinant not zero there; also Z m (x) is composed of elements 
analytic at x = a m or with a pole at that point. 

§ 11. Final Form of Solution. 

There is a certain lack of symmetry between the r6Ie of a h . . . , a m 
in the solution of the Riemann problem obtained in § 9, provided the 
given Cauchy matrices Z x (x — a^),. . ., Z m (x — a m ) were taken in the 
most general possible form. We shall now proceed to show that if 
Z m {x — o m ) be a properly chosen Cauchy matrix associated with the 
transformation T m , the equivalence of Y (x) and Z m (x) can be made 
proper at x = a m also. 

A form of reduction that is well-known suffices to establish this 
fact. 12 Nevertheless, inasmuch as a similar type of reduction is neces- 
sary later in the present paper, I give this reduction herewith. 

Let us begin with the matrix Y(x), obtained in § 9, which we may 
assume to be improperly equivalent at a m to the Cauchy matrix 
Zm (% — a m ). Now we have (§ 8) 

Z m (x — a m ) = I' (x — a m ) C 
so that from (30) 

Y(x)C-^ = Z m (x)I'(x-a m ), 

12 Plemelj, loc. cit., pp. 237-240. 
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The matrix of elements on the right can be written in the form 

(x— a m ) kl (an+b n (x— a m ) + ...),..., (x—a m ) k >>(a ln + 6 lB (x— a m ) + . . .) 

(33) 

(x—a m ) kl (a nl +b n i(x—a m ) + . ..),..., {x—a^) kn {a„ n +b m {x — a m ) + ■ ..) 

Here it is supposed that the exponents of the Cauchy matrix do not 
differ by integers; but a similar form can be found in all cases. In 
each column of Y(x) the highest possible power of a: — a m is exhibited 
which leaves the coefficients of this power analytic in character at 

X "m* 

If | ciij | 5^ we can write this matrix in the form A(x)I' (x — a m ), 
where I' (x) is the matrix (x k idij), and where A (x) is the matrix ob- 
tained from (33) by striking out the exhibited powers of a; — a. From 
the equation 

Y(x) = A{x)T{x-a m )C- x 

we see that Y(x) is properly equivalent to a Cauchy matrix belonging 
to T at x = a m . 

On the other hand if [ a^ | = we proceed as follows : It is readily 
verified that if Y(x) be one matrix solution of the Riemann prob- 
lem satisfying the relation of equivalence given in § 10, then 
DY (x) [ | D | ¥" 0] is also a solution. Consequently it is no restric- 
tion in the consideration of (33) to assume that 

«11 == «12- • • = llti = 0, 

since when | ay | = this relation may always be made to hold 
by multiplying on the left by a suitable matrix D. Denote by 
yiix),. . . ,y n (x) the elements of the first row after a factor (* — a m ) 1 
has been removed, where I is the exponent of the highest power of 
x — a m that may be taken out and leave the elements y\ (x), . . ., y n {x) 
of the respective forms 

(x— a m ) k ' (/i+ gi (x — a m ) + ...),. . ., (x — « m )M/n+ 9n{x — a m ) + . . .) 

It follows that /i, . . . , /„ are not all zero. 

We may now add constant multiples of y\ (x), . . ., y n (x) to the succes- 
sive rows of Y(x) and obtain a new matrix Y(x) . In fact this is equiva- 
lent to multiplying Y(x) on the left by a matrix 



P(x) = 
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1, 0, 

C2 , i, o 



(x — a m ) 1 



0, 1 



(a; — a m ) 1 



Inasmuch as | P(x) \ = 1, and the elements of P(x) are analytic save 
at x = a m , this modification cannot affect any of the properties already 
secured for Y(x). However by choosing c 2 , . . ., c m properly we can 
clearly make all the coefficients a 2 ,-, . . . a ni vanish at x = a m pro- 
vided fi ^ 0. In this manner we can increase an exponent k t by 1 
without altering the determinant of Y(x) except by a constant factor. 

Inasmuch as | Y(x) | does not vanish identically, a succession of 
steps of this type will finally bring to light a solution Y {x) for which 
Y(x) C has the form (33) and in addition | a# | y^ 0. When this 
stage is reached Y(x) will be properly equivalent to a Cauchy matrix 
belonging to T m at the point a m . 

When Y(x) has thus been given a normal form, it is the solution 
of a linear differential system (24) with regular singular points at 
x = a,i,...,a m and having no other singular points, as may be at once 
proved. 13 The elements of R(x) therefore have the form of rational 
functions whose numerators are polynomials in x of degree at most m — 2, 
and whose denominators are the product of (x — Oi), . . ., (x — a m ). 

§ 12. Irregular Singular Points and Canonical Systems. 

The Cauchy matrix is the simplest possible matrix of functions to 
which a matrix solution of a given linear differential system is properly 
equivalent at a regular singular point. Let us determine the simplest 
possible matrix Z(x) to which the matrix solution Y(x) of a given linear 
differential system (24), in which the elements of R(x) need not be 
rational, is properly equivalent at a prescribed irregular singular point. 
It is convenient to take this point to lie at infinity. If the highest 
order of any pole of an element of xR(x) at x = oo is p + 1 (p ^. 0), 
then p + 1 is said to be the rank of the singular point x = <*>. 

According to the results of § 10 we can find a matrix Z(x) which at 
x = co is properly equivalent to Y (x) and at another point x = 



13 Schlesinger, loc. cit. pp. 215-221. 
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is improperly equivalent to a Cauchy matrix which at x = undergoes 
a transformation inverse to that which Y (x) undergoes at x = oo . 
Here we take m = 2, a x = 0, a? = oo. The condition T2T1 = I is 
satisfied. 

By means of a modification precisely like that of § 11 we can make 
Z(x) properly equivalent to a suitable Cauchy matrix at x — 0, and 
yet preserve the other properties listed in § 10. 

Now consider 

ft (*)-^ **(*). 

Since | Z \ (x) | 5* for x ^ 0, oo , and since Z(x) and dZ(x)/dx 
undergo the same substitution about x = 0, the elements of ft(.r) are 
single-valued, and analytic for x ^ 0, oo . Since Z{x) is properly 
equivalent to a Cauchy matrix at x = 0, the elements of Ri(x) have 
poles of at most the first order at x = 0. 14 Moreover since Z{x) is 
properly equivalent to Y(x) at x = oo we have 

Z(x) = A(x) Y(x), 

where the elements of A (x) are analytic at x = oo and also | A(x) | ?* 
at x = oo . Therefore we obtain 



dZ(x)_ A{x) dYM + dA Y{x) = 



dx dx dx 

and 

. . dZ {x) „_, , , 
ftO) = —fa~ Z W 



^ (x) fi (a:) + ~ 
dx 



Y(x), 



~ A{x)R{x) + dJM 



A- 1 (x). 



Hence the elements of Ri(x) are analytic or have a pole of order not 
greater than p at x = oo . 

From this analysis it follows that xR\(x) is a matrix of polynomials 
of degree at most p + 1, so that Z(x) is itself the solution of a linear 
differential system 

(34) x~ = P(x)Z, 

where P(x) is a matrix of polynomials of degree at most p + 1. This 
is the canonical form of equation with an irregular singular point of 
rank p + 1 at x = oo . 

At a finite singular point x = a, the canonical system is of a type 
obtained from (34) by a transformation x' — a = l/x. 

14 The detailed proof is entirely similar to that given herewith to determine 
the nature of the elements of R (x) at x = oo . Cf. Schlesinger, loc. cit., 
pp. 143-144. 
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In order not to introduce artificial difficulties we shall consider the 
regular singular point to be of rank zero. To this case the above 
argument applies also, and the canonical system is that satisfied by a 
Cauchy matrix. 

If Y(x) is the matrix solution of a differential system (24) in which 
the elements of R(x) are rational, the above argument shows that at 
each of its singular points a h ... , a m , Y(z) is properly equivalent to 
the matrix solutions Z\(x), . . . , Z m (x — a m ) of canonical differential 
systems. 

Conversely the results of § 10 lead to the conclusion that given 
Ti,...,T m such that T m T m -i . . . Ti = I, and canonical differential sys- 
tems belonging to the singular points ai,...,a m with matrix solutions 
Z\ (x), . . ., Z m (x) undergoing a transformation to Z\ (x) TV, ... , Z m (x) T m 
at these points respectively, there will exist a matrix solution Y(x) 
of a rational differential system (24) which undergoes a transformation 
to Y{x) Ti as x makes a positive circuit of a { and which is properly 
equivalent to Z^x) at this point, for i = 1, . . ., m — 1, and properly or 
improperly equivalent to Z m (x) at x = a m . 

It is therefore essential to obtain a characterization of the matrix 
solution of a canonical system (34), and further to solve the associated 
inverse problem, before solving the general problem of characteriza- 
tion for a system (24) with irregular singular points. 

§ 13. The Problem of the Irregular Singular Point. 

In my paper referred to 15 I characterized the nature of the matrix 
solution of a canonical linear differential system (34), at least in the 
case that the roots of a certain characteristic equation were distinct; 
the case of equal roots introduces complications of an algebraical 
nature, and is put to one side in the present paper. 

The results which I obtained may be recapitulated as follows: 
If the singular point is taken at x = oo , there exists a formal matrix 
solution of (34) 

S(x) = (e v J^x k JSij(x)), 

x p+i x p 

(35) pj(x) = aj— — - + &•—+ • • • + hx, U = 1, • • -, n), 

su(x) = sn + s;/ 1 )- + . . ., (I, j = 1, . . ., n), 

is I, §§ 6, 7. 
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where | s^ | 5^ 0. The quantities <xi, . . ., <x„ are the roots of the char- 
acteristic equation alluded to, and the series s t j (x) are in general di- 
vergent. We shall assume for the time being that no three of the 
points ai, . . . , a n lie on the same straight line in the complex plane. 
Let now n, . . .,t n denote the N = n(n — 1) (p + 1) arguments 
in order of increasing angular magnitude such that for some j and 
k (j * k) 

(36) 9t { (<y — a k ) x p+I } = 0, arg x = r m . 

Here "9t" denotes "the real part of". Let us write tjv+i = t\ + 2tt, 
and let j m and k m denote the value of j and k corresponding to m, so 
ordered that the real part (36) changes from positive to negative 
as arg x increases through arg x = r m . There exist then N matrix 
solutions Zi(x), . . ., Z N (x) such that f or i = 1, . . . , N 

(37) Z m (x) ~ S (x), T m ^ arg x < T m+ i, 16 

and such that along arg a: = r m+l , Z m+i (x) and Z m (x) differ only in 
their j m th column, the j m th column of Z m+ i (x) being obtained from 
that of Z m (x) by the addition of the k m th column of Z m (x), affected 
with a suitable constant multiplier c m , to the j m th column. As a 
matter of definition we take 

(38) Z N+1 (x) = Z 1 (x)I', V = (e^-Hn). 

The proof of the existence of Z\{x), . . ., Z N (x) having these proper- 
ties can be directly based on the existence of a solution Z (x) asymp- 
totically represented by S (x) along every particular ray. 17 

The properties so far stated are characteristic of the behavior of 
the matrix solution not alone of a canonical system but of any system 
with singular point of rank p at x = 00 in the neighborhood 
of the singular point, 18 and are invariant under a transformation 
Y(x) = A (x) Yix) where A (x) is a matrix of elements analytic at the 
singular point in question and of determinant not zero there. 

When Z\ (x), . . . , 2jv+i (x) are in addition the solution of a canonical 
system, the matrices Z m (x) are analytic in the finite plane for i^O 



16 The relation "z (z)~s (x), arg x = a," means in the present paper that 
2 (x) is asymptotically represented by s (x) in some sector (however small) that 
includes arg x = a as an interior ray. This slight modification of the conven- 
tional meaning of the symbol " « " and its natural extension to matrices is 
convenient for the present paper. 

17 I, § 6. 

18 I, § 6. 
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and of determinant not zero. At x = these matrices are properly 
equivalent to a Cauchy matrix. These further facts come at once 
from the form of (34). 

Conversely let us prove that if Zj (x), . . ., Zjf (x) exist possessing 
such properties, they are all matrix solutions of one and the same 
canonical system (34). 19 

In the first place it is clear that the matrix 

P(x) = x ^^ Z m ~i (x) (m=l,...,N) 

is denned in each sector r m ^ arg x < r m+ i (i = 1, . . . , N) and hence 
is defined in the entire plane. Now Z m+ i (x) is obtained from Z m (x) 
by multiplication on the right by the matrix 

where C m is a matrix of zero elements except for the single element c m 
in the k m ih row and j m ih column. It appears therefore that P (x) is 
single- valued and analytic along each ray arg x = r m ; for we have 
along this ray 

^^W 1 (x) = dZ ^L[I + Cm ] [I + OJ-'Z m -i (x) 

- dx t* w- 

Moreover, on account of the equivalence of Z\(x), . . . ,Z m (x) to a Cauchy 
matrix at x = 0, the elements of xP (x) are analytic at x = 0. 

Secondly, in the neighborhood of x = <», the matrix P{x) is asymp- 
totic to 

(39) x — ^ — S~ l (x), Tm< arg x < r m+1 , 

where the meaning of the notation is manifest. It is to be recalled 
that the relation (37) holds (by definition) in a small sector including 
arg x = r m as an interior ray. This enables us to write 



dZ m (x) dS (x) . 

~dx~^-dT' r m ^ g x< Tm+ , 



20 



19 I stated this fact without proof earlier; see I, p. 468. 

20 Ford, loo. cit. 
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Thus P(x) is represented asymptotically in the complete vicinity of 
x = oo by (39). But we have formally 

(40) ^ = («,«{(«£) + 5>,,(,, + *|M}). 

If (35) and (40) be used in evaluating the expression (39), it is seen 
that each element of P(x) is given asymptotically by a power series 
in descending integral powers of x, with leading term in x p+1 or lower 
power of x. It follows that the elements of P{x) are analytic or have a 
pole of at most the (p+l)th order at x = <x> . Hence the elements of 
P{x) are polynomials of degree at most p + 1. 

The central problem of the irregular singular point is, for a given 
choice of pi(x), . . ., p n {x), r lt ...,r n and of c\, . . ., c# above described, 
to construct a matrix Y(x) with the above specified properties. 21 



§ 14. Solution of the Problem of § 13. 

In order to solve the problem just stated we make an application of 
the preliminary theorem of Part I, taking r = N/2 and for the curves 
Ci, . . . , Cjiv not the N/2 straight lines formed by the N rays arg x — r m 
(m = 1, . . . , N) but by the N rays arg x = r' m (m = 1, . . . , A T ), which 
are obtained from them by a slight rotation e in a clockwise direction 
in the complex x-plane. It is evident that if e be taken small 
enough we shall have 

Z m (x) eo S (x), r'm ^ arg x ^ r' m+ i, 

for m = 1, . . ., N. Furthermore we shall have for any m 

9* ( a im — «m) * p+1 > 0, arg X = r'm, 

for m = 1, . . . , N. This fact is essential to the solution. 

The matrices A\{x), . . . , A^{x) which are to be used in the applica- 
tion of the theorem are defined as follows. Write 

T{x) = (x p jWx r i6ij), 

and then put 

(41) A m {x) = T{x){I+ C m }T-> (x) (m = 1, . . ., N), 

21 Cf. I, § 7. 
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where x in A m (x) is taken along arg x = r' m , and the determinations 
of T(x) chosen are obtained from one another by allowing arg x to 
increase from n to t n . The matrix A m (x) thus defined is analytic 
along its line save at x = 0. Also A m (x) — J is a matrix whose only 
non-zero element is 

C m e p k m ( x )-PimWx r km- r im. 

Now we have 



v 



Pk m (?) — Vim (x) = (ak m — aj„) "j-r + (ft m — ftj— 

V + i ? 

+ • • • + (Xi m X; m )z, 

which quantity by definition of j m , k m has a negative real part for a; 
along the ray arg x = r' m , at least when | x | is sufficiently large. In 
consequence of the form of this non-zero element it is certain that 
A m (x) is unlimitedly differentiate along the ray at x = oo and be- 
haves there like the matrix /. As this is true along each of the rays, 
the matrices A m (x) clearly satisfy the conditions of the theorem in the 
vicinity of x — oo. Choose Ai(x), . . ., A$ N (x) along_the straight lines 
Ci, . . . , Ctf/2 as equal to the corresponding matrix A m (x) along either 
of its component rays outside of some circle | x \ = r, within which 
they are chosen so as to satisfy the conditions of the preliminary 
theorem. Since A 1 (x), . . . , A N (x) are of the nature above described 
at x = oo , the matrices A x (x), . . . ,A N/2 (x) are unlimitedly differen- 
tiate, and satisfy the permutability condition (see (21)) of the pre- 
liminary theorem. It is therefore possible to make such a choice 
of A x (x), . . . A N/i (x). 22 

It follows by the theorem that there exists a matrix $(cc) of determi- 
nant nowhere zero save possibly at x = a = 0, analytic except along 
the rays arg x = rj and such that 

(42) x ^ m+ $ ( X ) = [Jta _ $ {x)] A m (x m ), 

where x m is a point of arg x = r m ' for which | x m \>. r. It also follows 
from the theorem that each element of <J>(a:) is represented asymptoti- 
cally by a series in negative powers of x in each sector (r' m , r' m+1 ) ; 
since A m (x) ~ /, these series are the same in all the sectors and we may 
write 

(43) *(*)~M*)), 

where s^x) is a series of negative powers of x in which the determinant 
of the constant terms is not zero. 

22 Cf. §9. 
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Now write for m = 1, . . . , A r + 1 

Z m (x) = $ (ar) T (x) (t'„ ^ arg z ^ r' m+ i). 
From (41) and (42) there results along arg x = r' m+ i 
Z w (x) = Z m (x) [I +C m ], | x | ^ r. 

The functions Z m (a;) may accordingly be continued analytically across 
each ray arg x = T m and represent matrices analytic for | x | ^ r, of 
determinant not zero. The relation (43) leads to the conclusion that 

Z m (x) ~ S(x) (r' m ^ arg a- ^T' m+ i)- 

where S(x) is of the desired form (35), and thence to the conclusion 
that this asymptotic representation is valid for r m ^ arg x < r m +\. 
The relation between Z^+i (x) and Z x (a-) is that stated in (38). 

In § 13 it was shown that, if matrices Zi(x), . . . , Zjf(x) had the above 
properties and the further properties that they were analytic in the 
finite plane of determinant not zero for x ?* 0, and at x = were 
properly equivalent to a Cauchy matrix, then these matrices were 
solutions of a canonical system (34) with irregular singular point 
at x = oo . The same arguments can be used to establish that the 
Z\ (x), . . .,Z N (x) before us are solutions of a differential system (24) 
having coefficients rational in character at x = oo , with poles of order 
not more than p. But we have proved in § 12 that the matrix solu- 
tion of such an equation is properly equivalent to that of a canonical 
linear differential system at x = oo . Consequently a transformation 
Z(x) = A (x)Z (x), where A(x) is a matrix of functions analytic in 
character at x = oo of determinant not zero there, makes Zix) the 
solution of such a canonical system. In particular the matrices 

\ (x) = A (x) Z x (a-), , Z N (x) = A (z) Z N (x) 

form the solution of our problem, a fact which is apparent if we 
note that S(x) = A(x) S(x) has the same form as S(x). 

It has therefore been completely established under the stated 
restrictions that the characteristic constants which occur in the 
characterization of the matrix solutions of a canonical linear differen- 
tial system can be chosen at pleasure. 

The restriction that no three of the quantities m, . . . , o„ shall lie 
on a straight line is not essential, for if it is not satisfied and if no two 
of the polynomials pi{x), . . ., p n (x) are identical it will be possible to 
replace the rays n, . . ., tn which have coalesced by an equal number 
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of curved rays so chosen that the real part of only one of the differences 
Pi (%) — Pj (•**) changes sign along the ray, and thence to apply the 
preliminary theorem in much the same fashion as before. If two or 
more polynomials pi (x), . . ., p n {x) are identical it is not necessary to 
modify the nature of the rays. 

Furthermore it would be possible to construct an analogous existence 
proof when S(x) has for its elements so-called anormal series. 

In other words, the results here obtained are of an entirely general 
nature. 

§ 15. The Generalized Riemann Problem. 

The problem which I proposed in my paper on singular points (I) 
was the following: " To construct a system of linear differential equations 
of the first order with prescribed singular points 

a'i, X2, . . . x m , x m+ i = co 
of respective ranks 

°u • • • > <7m+l, 

and with a given monodromic group, the characteristic constants being 
assigned for each singular point." This problem is virtually solved by 
what precedes. 

It is of course understood that certain obvious conditions of com- 
patibility are satisfied, the first being the one already noted, 

-t m+l-l m ... 1\ — 1. 

However a second necessary condition must also be imposed. Take 
any assigned point x = x { . The solution Z x (x) appertaining to this 
point (see § 14) is transformed successively into 

Z 2 (x) (1 + CO" 1 , Z z (x) (I + C 2 )-' (/ + d)-', 

...,z N+1 (x)(i+c N yi... {i + c 1 y\ 

as x passes over the rays r x , . ■ . r N respectively. Hence after a com- 
plete circuit of a,-, Z t (x) alters to 

' Z 1 (x)I'(I+C N )-K.. (J+d), 

where the matrix T t of transformation is explicitly determined in 
terms of the characteristic constants. But in order for this set of 
characteristic constants to be possible, some solution Z (x) = Z x (x) C 
must undergo precisely the transformation by T h i. e., 

Ti^CTiC-* (i=l,...,m). 
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This second condition is satisfied if the elementary divisors associated 
with Ti and T t are the same. 

With the understanding that the two conditions of compatibility are 
satisfied, we can assert that there exist m + 1 canonical systems with 
the prescribed characteristic constants at X\, . . ., x m+1 and with solu- 
tions Vi (x), . . ., V m+1 (x) undergoing a transformation to V\ (x) T\,..., 
V m +\T m +\ as x makes a positive circuit of x\, . . . x m+ i respectively. 
Hence there will exist a matrix Y (x) of elements analytic and of de- 
terminant not zero save at Xi, . . . , x m+ i, properly equivalent to V t (x) at 
x = a,- (i = 1, . . . ,m), and properly or improperly equivalent to V m+ i (x) 
at x m+ i = oo . This follows upon application of the results of § 10. 
The matrix Y (x) thus obtained is a matrix solution of a differential 
system of the form required and has the monodromic group and 
characteristic constants required, if the equivalence at x m+ i = » be 
proper. This is an immediate consequence of the results of §§ 13, 14. 

If this is not the case it is easy to show that if one merely increases 
the characteristic constants n, . . . ,r n for the singular point x m+i = oo 
by suitable integers this equivalence at infinity becomes proper. For 
consider the matrix S(x) corresponding to the matrix solution of the 
equivalent canonical system at x = c© . The corresponding matrix 
S(x) for Y{x) is of the form A{x) S(x) where the elements of A(x) are 
rational in character at x = °° . Hence if a proper choice of C be 
made, the matrix S (x) C may be written (compare (33), § 11). 



gPi( x )X r i 


(- + £+■■■} 


. . ., e p n (x ^X rn 


\ X J 


e v ^x f i 


(- + £+•■■) 


. . ., e p n^x r » 


(a m + b f+...) 



where n, . . . , r n differ from r\, . . . , r n by integers. By a process of re- 
duction exactly like that given in § 11 this matrix may be replaced by 
a series of matrices of the same form in which the exponents r h . . . , r n 
are increased so long as | ay | 5^ 0. We conceive of Y (x) as affected 
by the same series of operations, which have no effect on the proper- 
ties above specified. Furthermore the transformed Y(x) and S(x) 
stand always in the same relation to each other after as before trans- 
formation. If the process comes to an end so that | a y - 1 9^ 0, the 
linear differential system with matrix solution Y(x) will have a singu- 
lar point of rank q m+ \ at x = 00 (compare § 13, pp. 547-548) and the 
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associated properly equivalent canonical system has characteristic 
constants only modified as stated. 
Now we have 

I Y(x) J ~ \S(x)\ 

in the complete vicinity of x = go , which is made up of the sectors 
(ti,t 2 ), . . . (tjv,tjv+i). It follows that the series for | S (x) | converges 
and hence must be of the form 



e Pl<*)+ . +P n (z) . X V L _|_ \. + _ . , \ a ^Q 
\ X 1 



so that always 

n+ . . . + r n ^p. 

But the above reductions increase r\ + ... + r n and do not affect 
| Six) | or the value of p, and so must terminate. 

We can state then that either a solution of the stated problem, or of 
a modified problem in which the constants r\, . . ., r n of one of the singular 
points are altered to n + l h . . . , r n + l n respectively, where I., . . .,l n are 
integers, will exist. 

The matrix Y{x) thus obtained is not always unique. The most 
general determination is however of the form P(x) Y(x) where P(x) 
is a matrix of polynomials of constant determinant which fulfills other 
conditions. Thus the notion of "primitive systems" admits of 
extension to the case of irregular singular points. 23 

Part III; The Linear Difference Equation Problem. 

§ 16. Formulation of the Problem. 
Let 

(44) Y(x+ 1) = Q(x)Y(x) 

be a linear difference system in which the elements of Q(x) are poly- 
nomials of degree ju in x. 2i In my earlier paper on linear difference 
equations I demonstrated that, at least if the above equation admits 
a formal matrix solution 

(45) S(x) = [x^{ P je-^ x x r 'sij{x)}, 

in which s tj (x) is a power series proceeding according to negative 
powers of x with the determinant of the leading coefficients not zero, 



23 Cf. Plemelj, loo. cit. pp. 240-245. Like theorems may be proved in a 
similar manner here. 

24 Essentially the most general linear difference system with rational coef- 
ficients may be reduced to this form; see II, § 5. 
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there exist two matrix solutions Y~ (x) and Y + (x), with elements 
analytic in the finite plane save for poles, such that Y~ (x) co S (x) in 
any left half plane and Y + (x) <» S (x) in any right half plane. 25 The 
existence of such a solution was proved by Norlund and Galbrun by 
methods based on the Laplace transformation somewhat earlier. 26 
These matrices Y~ (x) and Y + (x) are connected by a relation 

(46) Y~ 0) = 7+ (.r) P (x) 

where P (x) is evidently a matrix of periodic functions of period 1. 

From the form of (44) it appears that Y~ (x) is a matrix of entire 
functions, while Y + (x) is analytic save for poles. 

In my paper I determined explicitly the nature of the elements 
Pij (x) of P (x) to be the following: 

(47) 
r Vii ( x ) = l -f afih 2 * ^' lx +...+ cue 2 *(K-» v -' x + e 2 "' ^e**» v ~' * 

(i= 1, . . ., n) 

P ij(x) = e^'i ^- lx [cijW + . . . + c»/*<"-« ^~ lx ] 

(i ^ j; i,j = 1, • • -, n), 
Here \j stands for the least integer as great as 

(48) 9t(^=Pogffi-lo gP ,])- 

An analogous determination in certain cases at about the same time 
was made by Norlund (loc. cit.). 

It is not difficult to show that, if Y~ (x) and Y + (x) have the properties 
above outlined, then conversely they are solutions of a linear difference 
system (44) in which the elements of Q (x) are rational if not polyno- 
mial. 27 For this reason the constants pj, Tj, c f / w may be called the 
characteristic constants of Y~ (x) and Y + (x). 

This characterization suggested to me the following problem: To 
construct a linear difference system (44) with assigned characteristic 
constants in which the elements of Q(x) are polynomials in x of 
degree not greater than ju. 

25 These matrices Y~(x) and Y+(x) correspond to 0(x) and H (x) of II. In 
certain cases it may be necessary to consider half planes not bounded by a 
vertical line. I refer to this possibility later. 

26 Norlund, Dissertation, Copenhagen (1911); Galbrun, Dissertation, Paris 
(1910). 

27 See II, § 7. 



BIEKHOFF. — THE GENERALIZED KIEMANN PEOBLEM. 555 



§ 17. Solution of the Problem of § 16. 

In order to treat the problem of § 16 we apply the preliminary 
theorem. We shall take r = 1, and take Ci to be the axis of imagi- 
naries in the complex plane unless | P(x) j = at a point of that axis. 

The matrix Ai (x) is taken equal to 

T 0) P 0) T- 1 (x), T (x) = {x" x (pjtr")* x r >5i,), 

except near to x = 0, where it is chosen in any way so as to satisfy the 
restrictions of the preliminary theorem there (compare § 9). Since 
the elements of T{x) are in general multiple-valued functions of x, 
it is necessary to specify which branch of T{x) to select. We shall 
choose a continuous branch of T(x) in the right half plane, and a 
continuous branch in the left half plane in such a way that these 
branches coincide along the upper half of the axis of imaginaries. 
The first factor T(x) in the expression for Ai(x), will be identified 
with the first of these branches, and the last factor T~ x {x) will be the 
inverse of the second of these branches. 

It is therefore clear that, along the upper half of the axis of imagina- 
ries, the element in the ith row and jth column of A\ (x) is, for i 9^ j, 

e 2*\a <-i*pX pf x x n-rj f^m + . . . + Cij ^-i) e 2^-i) V=I*j 

while the diagonal elements are the same as for P(x). But by defini- 
tion Of Xy, 

(49) 1 > 9t (* {j - 7~j=^ (log w - logPi)) ^ 0. 

Let us exclude at present the case of the equality sign; the element of 
Ai(x) in the ith row and jth column (i ^ j) will therefore vanish to 
infinite order together with its derivatives as x goes to infinity along 
the upper half of the axis of imaginaries. The diagonal elements 
diminished by 1 have the same properties. 

Hence we have Ai(x) «> I along the upper half of the axis of imagi- 
naries, while all the derivative matrices of Ai(x) tend to matrices of 
zero elements as x becomes infinite. If Ai(x) has this character along 
the lower half of the axis also, it is clear that this matrix satisfies all 
the restrictions imposed in the theorem. 

Let us demonstrate that such is actually the case. The determina- 
tion of T(x) on the left-hand side of the lower half of the axis of imagi- 



556 PROCEEDINGS OF THE AMERICAN ACADEMY. 

naries is obtained from that on the right-hand side by a complete 
positive circuit of x = 0, during which T(x) changes from 

(x^ (pj<r") x x r i5ij) to (e 2ir * VIi V* V=I r ix^ (pje~") x x r i5y) . 

The ith. diagonal element of Ai(x) may now be written 

Ci .(0) e -27r iJ-lr Je -2TV V=> x _]____)_ ^ 

while the element in the t'th row and jth column of A (x) (i 9^ f) may 
be written 

Bearing (49) in mind we readily perceive that Ai(x) does have the 
indicated properties along the lower half of the axis of imaginaries. 

According to the preliminary theorem we can then determine a 
matrix $(a;) such that 

(50) ,2^*(*) = [A*(*)Mi(z), 

where xi is a point of the axis of imaginaries and the approach is from 
the left-hand and right-hand side of that axis respectively. If we take 
x = a = 0, the determinant of $(x) is not zero in the finite plane except 
at x = possibly, and the elements of this matrix are analytic at any 
point not on the axis. Along the axis as defined from either side these 
elements have continuous derivatives of all order, and will be analytic 
at more than certain distance d from the origin. In the vicinity of 
x = 00 , $(.*•) is represented asymptotically by a matrix of series in l/.r 
with determinant of leading coefficients not zero. This matrix is the 
same on either side of the axis, since A (x) «> J. 

Let us denote $(x) by U + (x) for x in the right half plane and by 
U~ (x) for x in the left half plane, and write 

y+ (a:) = [7+ (x) T (x), Y~ (x) = U~ (x) T (x). 
From equation (50) we see then that 

(51) Y-(x) = ?+(x)P(x) 

for I x I > r along the axis of imaginaries. From the asymptotic 
form of U+(x) and U~{x) at x = <x> we obtain 

Y-(x) ~ S(x), F + (x) ~ S(x), 

in the left and right half plane, where S(x) is of the same form as S(x) 
above. The relation (51) shows that Y~(x) is composed of elements, 
analytic in the right half plane. 
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Let us now apply the preliminary theorem a second time, taking 
r = 1, and for C\ a circle with center at the origin and radius so large 
as to include within it all those points of the axis of imaginaries at 
which an element of Ai(x) as chosen above is not analytic, and also 
so as not to pass through a zero of [ Y~(x) \. 

In this second application of the theorem we choose Ai(x) to be 

[F"(.x-)] _1 , and in this way satisfy the restrictions of the theorem. 
Furthermore let us take a = oo . 

Along the circle Ci we have for the solution $(«) 

(52) ,2S + * (*) = [,!£- $ (*)] [Y- (x 1 )}-\ 

where the approach to the point x\ of C is from without and within C 
respectively. Now write 

Y- (x) = $ 0) Y- (x), Y+ (x) = $ (a;) 7+ (a;) 

for x outside of C\. It follows that Y~(x) is composed of elements 
analytic in this region; also along C\, Y'(x) coincides with the inner 
determination of $(x) by (52), so that the elements of Y~(x) are also 
analytic within and on C\. Hence Y~(x) is a matrix of entire functions. 
Similar considerations show that the elements of Y + (x) are analytic 
in the right half plane like the elements of Y + (x). 

At x = oo , Y~(x) and Y + (x) are asymptotically represented by a 
matrix S(x) in which however r h . . ., r n are not necessarily the same 
as in S(x), and in which the determinant of the leading coefficients 
may be zero. This results from the fact that the elements of $(x) 
are rational in character at x = as and in consequence can be ex- 
panded in convergent series in descending integral powers of x. In- 
asmuch as we have | S (x) | = | $ (x) | • | S (a;) (, it is also true that 
| S (x) | cannot reduce formally to zero. 

Finally from (51) and (52) we infer that 

(53) F-(a-) = Y+(x)P(x). 

A first conclusion to be derived is that Y~ (x) •*> S (x) in any left half 
plane, and that also Y + (x) «> S (x) in any right half plane. In fact 
we have already determined the asymptotic form of P (a;), and this 
known form combined with the known asymptotic form of Y + (x) in 
the right half plane gives us the form of Y~(x) in the part of the plane 
to the left of any line parallel to the axis of imaginaries; a similar 
remark applies to the asymptotic form of F + (a-) in any right half 
plane. 
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One further remark comes in appropriately at this point. The 
definition of Y (x) and Y~(x) ensured that | Y~ (x) \ and [ Y + (x) [ do 
not vanish to the left or right of the axis of imaginaries respectively, 
save possibly for x = oo . Hence | Y~(x) | and | Y + {x) | do not vanish 
in these left or right half planes respectively. 

We may now enter upon a series of modifications of Y(x) which will 
preserve the above stated properties and secure in addition that the 
determinant of the leading coefficients of S(x) is not zero. To this 
end we write the matrix S(x) in the form 

x»* ( Pl e-») x x r < (a n + ^ + . . X . . ., x" x ( Pn e-») x x"n L ln +~+ ■■ 



Onl , \ „ T / ,.w - ( , b n 



xi* ( Pl e-") x x r H a„i+ —+...,.. ., x^ { Pn <r») x x r » a nn +-?! + . 



x 



and carry out reductions parallel to those given in § 11. The same 
reductions are supposed to be simultaneously effected upon Y~(x) and 
Y + (x) (compare § 15). This set of reductions will terminate, since 
| S(x) 1 does not vanish identically, and, when it does, F~(.x) and Y + (x) 
will have the desired additional property. 
Consider now the matrix 

Q (x) =Y-(x+ 1) [Y~ (x)]- 1 = 7+ (a: + 1) [F+ (*)]->. 

These two forms for Q(x) are equal in virtue of (53) and yield us at 
once the asymptotic form of the elements of Q(x) in the complete 
vicinity of x = °o as descending power series in x with leading term of 
the /uth degree at most. 28 Hence the elements of Q (x) are rational 
in character at a: = <x> , with a pole of at most order /x there. 

In the finite plane to the left or along the axis of imaginaries the 
first expression for Q(x) shows that the elements of Q(x) are analytic 
without exception. It will be recalled that Y~(x) is a matrix of entire 
function of determinant not zero in the left half plane inclusive of the 
axis of imaginaries. On the other hand the elements of Q(x) are ana- 
lytic to the right of the axis of imaginaries, as the second form shows. 

Accordingly Q (x) is a matrix of polynomials of degree at most ju 
and Y~(x), Y + (x) are solutions of the rational difference system (44). 

Finally it may be observed that in case | P(x) | = along the axis 
of imaginaries, a parallel line may be used to take the same role; or 

28 Cf. I, § 7. 
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indeed any simple analytic curve without a horizontal tangent and 
with vertical asymptote, provided that | P(x) | ^ along the curve. 

If the equality sign obtains in (49) it will be necessary to employ 
a curve with asymptote not quite in the vertical and to employ half 
planes not bounded by a vertical line. 

It is also possible to replace S(x) by certain anormal forms, 29 and 
thus extend the above results to the most general case. 

Our conclusion may be summed up as follows: There exists a linear 
difference system (44) with matrix solutions Y~(x), Y + (x) which either 
possesses prescribed characteristic constants pj,rj, c t f k \ or else constants 
Pi' r i "f" h' c v k) where l\, . . .,l n are integers. For an arbitrary curve 
which meets each line parallel to the real axis only once, having a ver- 
tical asymptote, and which does not pass through a point \ P(x) [ = 0, 
there exist such matrices Y~(x), Y + (x) with the further property that 
| Y~(x) \=A to the left of the curve while the elements of Y + (x) are ana- 
lytic and | Y + (x) | ^ to the right of the curve. 

It is worthy of note that this last property determines the location 
of the poles of the elements of Y + (x) completely: namely, they occur 
to the left of the curve and at the points for which | P(x) | = 0. This 
appears from the formula 

F+(x) = Y-(x)P- J (x), 

which also permits us to affirm that the precise maximum order of pole 
of any element of Y + (x) is the order of the zero of | P(x) |. 

Part IV: The Linear ^-Difference Equation Problem. 

§ 18. On Linear q-Difference Equations. 

A linear g-difference system may be written 

(54) Y (qx) = Q(x)Y (x) [ q [ > 1, 

where Q(x) is a matrix of polynomials of degree n or less, in analogy 
with the normal form (46) of linear difference systems. The apparently 
more general case in which the elements of Q(x) are rational in x may 
be reduced to this form readily. Let the least common denominator 
of the elements of Q(x), written as quotients of relatively prime poly- 
nomials, be 

(x — ai) . . . (x — ai) 

29 Analogous to the anormal series for linear differential equations. These 
forms have recently been obtained by Mr. P. M. Batchelder. 
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and let g t (x) be a solution of the g-difference equation of the type (54) 

(55) g (qx) = (x - m) g (x) 
for m = a,-. If one takes for new variable 

Y(x) = g 1 (x)...g l (x)Y(x), 

a new matrix equation (54) in Y(x) is obtained with Q(x) polynomial 
in x. 

Let us write 

= log a; 
log?' 

In terms of this new variable a solution of (55) for m = is 

qi(f—t)_ 

For ra^O, the transformation 

a; = war, y (x) = e" v lo 8 9 miog « y (a;) 

takes (55) to the normal form 

(56) y(qx) = (1 — x)y(x). 
Two solutions of this equation are 



(57) 



as one may verify by direct substitution. The function y x (a-) plays 
the same role for the linear ^-difference equations as the gamma func- 
tion does in the theory of linear difference equations. I have men- 
tioned these functions in order to supply an example later. 

The fundamental existence theorems for linear g-difference equa- 
tions are essentially a consequence of the work of Grevy 30 and Leau. 31 
The first complete treatment has been given by Carmichael, 32 and the 

30 Paris thesis, 1894. 31 Paris thesis, 1897. 

32 Am. Jour. Math., 34, 147-168 (1912). 
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result may be expressed as follows : There exist in general two matrix 
solutions 

{ F (x) =(x»iaij(x)) 
(58) „ ( 

where each function ay (x) is analytic at x = 0, and each function by (x) 
is analytic at x = oo ; and where furthermore the determinants of the 
leading coefficients of a { j (x) and 6 i; - (x) at x = and x = oo respec- 
tively are not zero. It is only the case when such series exist that 
will be here considered. 

It follows at once from (54) that Fo (a;) is a matrix of functions ana- 
lytic for x 5^ 0, oo , and that Y~ (a-) is a matrix of functions analytic in 
the finite plane except for poles when x 7^ 0. Further, if we write 

F (x) = Y x (x)P(x), 

then P(x) is a matrix of functions analytic for x ^ 0, oo , and possessing 
the property that P (qx) = P (x). 33 These properties are in close 
analogy with the properties for a linear difference system, to which 
indeed (54) reduces formally by the substitution x = q l . 

I propose now to determine completely the nature of P(x), as I 
have done for the analogous functions P(x) associated with the linear 
difference system ; it is the doubly periodic functions which enter here 
instead of the simply periodic functions. By means of this determina- 
tion it will be possible for us to state the problem which, for this field, 
is analogous to the problems above treated for linear differential and 
linear difference equations. 

§ 19. On the Matrix P(x). 
Let us make the transformation x = q l and write 
P(x) = P(t). 

The function Pit) is a single-valued function of t analytic save for 
poles; for, this transformation takes the Riemann surface of infinitely 
many leaves, with logarithmic branch points at x = and x = oo in 
a one-to-one and conformal manner into the £-plane. 

Let us conceive of the £-plane as divided into parallelograms which 

33 Cf . Carmichael, loc. cit., p. 159. 



562 PROCEEDINGS OF THE AMERICAN ACADEMY. 

belong to the periods co = 1, a/ = 2W — l/logg, and let ABCD 
be a parallelogram with vertices 

_i_ i , 1 , 2xV^ l 2xV^l 

Zfl, «o + 1, X + 1 H j , .To -\ -. , 

log? _ logs 

respectively. At homologous points of BC and AD, P(t) has the same 
value, since P(t + 1) =P(t). 

To obtain the relation between P(t) at homologous points of AB 
and DC we consider first the matrix 

(59) P(x)=Y x -i(x)Y (x). 

It is apparent from the form of the elements of Y (x) near x = as 
given by (58) that, if a positive circuit of x = be made, Y (x) 
will change to Yo(x) K, where 

likewise, upon a similar circuit, Y x (x) will change to 
where 

If these modified matrices be substituted in (59) we obtain the form 
which P(x) assumes after x has made a positive circuit of the origin. 
This is 

(_ lye- 2 ™ V=I 'J^Z-'P Or) K. 

If therefore pij (x) denotes the element in the ith row and jth 
column of P{x), such a circuit modifies p t j (x) to 

(_ iy e ~ 2 ™ V3i 'J^ie^i+tf ttpq (x). 

But this circuit in the a>plane corresponds to a passage from a point 
of AB in the i-plane to the homologous point of DC; in this way, let- 
ting Pij{t) stand for the element of P(t) in the ith row and jth column, 
we find 

(60) fait + 2 - J ^~) = (- D* 6 " 2 " ^ '^^'W <~ x m 



(63) \ 
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We have also seen that 

(61) pij(t+l) =pij(t). 

Now let us attempt to satisfy (60) and (61) by writing 

(62) pij (0 = ce a ' !+M <r (t — oi) . . . <r (t — a„), 

where a{t) is the Weierstrass sigma function belonging to the periods 

CO = 1, CO = — -. , 

logg 
which satisfies the relations 

„ (< + ^) = _/(^), (0 , 

log? J 

<r (t + 1) = - e"«+» <r (0, r? 2 ^ ~ 1 - V = 2tt V~l. 

logg 

The above choice of co and co' meets the requirement 91 (co'/co) > 0, 
since | q | > 1. 

A direct substitution into (61) determines 

m 

(64) a=— r,%, b = ^Eai+juW— l + 2farV— 1, 

in which k denotes any integer. If these values of a and b are taken 
and a direct substitution is made in (60), there is obtained the further 
condition 

(65) 2>H ; = a + pj+ I , , 

,ti log? logg 

in which I denotes an integer. These conditions are necessary and 
sufficient that an expression of the form (62) shall satisfy both (60) 

and (61). If the value for £ a > deduced from (65) be used in (64) 

the expression for b simplifies to 

V (o-j + pj + V) — r{ ( - + k). 

But if one adds or subtracts a period to a t the precise effect is to alter 
k or I by an integer. It is therefore always permissible to take k = I = 
and to write 

— w - , r , ., ... iVi 

(66) p««) = cue 2 + L ,l< " +P -' J 2 J< a (f _ ai (v-,) ...„(*_ 0j>( y)) f 
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provided that 

» U.T V 1 

(67) Z <*x (i - y) = °i + Pi i 

xti log? 

This is of course under the assumption that Pij(t) may be represented 
in the form (62). But this fact may be proved at once. For let \f/{t) 
be any function which satisfies (60) and (61), and <f>(t) the particular 
one above obtained. The function \p (t)/<j> (t) is doubly periodic, 
analytic save for poles, and can therefore be expressed as a quotient 
of products of sigma functions 

r (r(t — yi).. .a(t—y k ) _ yR 

But this quotient when multiplied by <j> (t), which is expressed as a 
product of sigma functions, must yield \p (t)> an entire function. This 
necessitates that for each zero of a (t — 7,-) in the numerator 
(i = 1, ...,/*) there must be a congruent zero x = fy of cr (x — /3 ; -) in 
the denominator. Such pairs of corresponding factors may be com- 
bined leaving only an exponential factor e ct+d . Thus \p (f) appears in 
the same form as <j> (i). 

Our result is therefore that the element p { j (x) of P(x) is of the form 
(66) (t = logo:/logg), where pi, . . .,p„, <ti, . . ., <r„ are the constants 
that appear in the series representations (58), and where the condi- 
tions (67) are fulfilled. 

It is interesting to apply these results to the equation (56) in which 
n = 1, n = 1, pi = 0, <ti = x V — l/log?, and Y (x) and Y x (x) reduce 
respectively to y (x) and y x (x) defined in (57). In this case we have 
therefore 

(68) 2/0 (x) = 2/oo (x) p (x) 
where 

The constant c may be determined by writing (68) in the form 

yo(aO = [(*-l) »«(*)] [^3] 

and allowing x to approach 1. Since o-(O) = 0,</(0) = 1, this gives 
(see (57)) 

(69) c(l-J) ! (l- ? „ 
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The relation (68) with the explicit values of y (x), y x (x) and p (x) 
substituted in is essentially one of the fundamental product formulas 
for the sigma function. 

§ 20. The q-Difference Equation Problem. 

It is now easy to show that conversely if To 0*0 > J» ( x ) are matrices 
of functions of the form (58) in the vicinity of x = and x = oo 
respectively, analytic for x 7^ 0, 00 , save for poles, and if the matrix 
P(x), defined by the relation Y (x) = Y x (x)P(x), is composed of 
elements Pij(x) which are left unchanged when x is replaced by qx, 
then Fo (x) and Y^ (x) are matrix solutions of a linear g-difference 
system (54) with rational coefficients. In fact, if we write 

Q (x) = F (qx) F„-' (*) = Y x (qx) Y^ (x), 

it is seen at once that for x t 6 0, 00 , the only singularities of Q(x) 
are poles, while the first and second of these forms for Q(x) ensure that 
Q(x) is composed of elements analytic at x = and with a pole of at 
most order n at x = 00 if not analytic there. This suffices to establish 
the fact that the elements of Q(x) are rational. 34 In order to conclude 
further that the elements of Q(x) are polynomials of degree n it is 
sufficient to know that the plane may be divided into two parts by a 
loop about x = meeting each equiangular spiral or radial line 

(70) 6 = c + j — t~ log r (r, 6, polar coordinates) 

only once and not passing through a point [ P(x) \ = 0, such that the 
elements of Y (x) are analytic and |F («)| is not zero within or along 
the loop, while the elements of Y x (x) are analytic and | F x (x) | is 
not zero outside the loop. Under these conditions the first expression 
for Q (x) makes it evident that its elements are analytic within or 
along the loop, and the second expression makes it clear that the same 
is true without the loop, since if a; is a point without the loop so is qx. 
It is natural to term the 2re constants pj, cxj and the n 2 (p. + 1) con- 
stants c t j, a\- l ' j \ . . ., a/ l ' j) the characteristic constants. These constants 
are not all independent, since there are n 2 relations between the con- 
stants ag/ 1 '^ of the type (67). Furthermore the constants c i; - are not 
uniquely determined by the given g-difference system. For, any ith 
column of Yq (x) is only determined up to a constant factor /,■ and like- 

34 Compare II, § 7. 
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wise any ith. column of Y x (x) is only determined up to a constant 
factor h ( . This fact appears from (58). The effect is to allow one 
to replace c i; - in Pijix) by ffij/hi, and thus vary at will 2» — 1 of the 
characteristic constants c i; -. There remain then essentially only 

2 n + n 2 (p. + 1) — n 2 — (2« — 1) or n 2 ju + 1 

characteristic constants. It is readily verified that these constants 
are all invariant under a transformation Y(x) = CY(x), where C is 
an arbitrary matrix of constants. 

But the equation (54) involves w 2 (/x +1) arbitrary coefficients in 
Q(x) of which there are n 2 (fj. + 1) — (n 2 — 1) or raV + 1 invariants 
under the same linear transformation. Hence we have found as many 
invariants for the ^-difference system as invariant characteristic 
numbers for the solutions. 

We are thus led to formulate the following problem: To construct 
a (/-difference system (54), with coefficients of degree not greater than 
ju in x, having any assigned set of characteristic constants. 

§ 21. Solution of the Problem of § 20. 

Here also we shall employ the preliminary theorem, but the appli- 
cation of it is even simpler than in the earlier cases. 

The conclusion that we shall derive is the following: There exists 
a linear q-difference system (54) with the matrix solutions Y$ (x), Y x (x) 
either possessing prescribed characteristic constants p i} cr,-, c t j, a-c"", . . . , a^ %,1) 
or else constants pj, c,- + l- p c i; -, a} Ul) , . . ., a^ ]) , where l u . . .,l n are integers. 
For an arbitrary loop about x = which cuts each spiral (70) only once 
and does not pass through a point | P(x) | = 0, there exist matrices 
Y (x), Y x (x) with the further property that \ Y (x) \ ^ within or along 
the loop while the elements of Y x (x) are analytic and \ Y x (x) | is not 
zero without the loop. 

Let Ci be a specified loop of this description which may be taken to 
be analytic. 36 We may take the matrix A\ (x) of the preliminary 
theorem to be 

T x (x)P(x)T -i(x), 

where ^«_ n 

T (x) = (x'idij), T x (x) = c^ (x-'iSa), 

since the elements of Ji(x) are single-valued and analytic along C t 
by (60). 

36 It is only a question as to how the loop weaves among the zeros of \P(x)\. 
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According to the theorem there exists a matrix <£(x) such that 

(71) ,!£+*(*)= [ x l ^ 1 ^(x))A 1 (x 1 ), 

where xi is an arbitrary point of C\ and the + and — signs denote 
approach to x\ from within and without C\ respectively; the matrix 
&(x) possesses certain other properties: it is composed of elements 
analytic for x not on Ci, save at a point a which we shall take to be at 
infinity; furthermore its determinant does not vanish in the finite 
plane. 

Let us denote by Uq (x) the matrix $ (x) within C h and its 
analytic extension across &; similarly by U x (x) let us denote the 
matrix $(x) without Ci and its analytic extension across C\. Consider 
then the matrices 

Y - (x) = U (x) To (x), Y x (x) = U x (x) T x (x). 

From (71) we obtain at once 

(72) Y -(x)=Y x (x)P(x), 

and prove without difficulty that Y ~(x) and F M +(x) as thus defined 
have the characteristics demanded save possibly that Y x (x) may not 
be precisely of the form (58) at x = oo , as it would be if $(x) were com- 
posed of elements analytic at x = oo and if also | $(.i-) | were different 
from zero at x = oo . Nevertheless one can always write F M (x) in 
the form 



1 (f-t) 



fflu + — + ...),.. . x-"n I a ln + 



X J \ X 

am +—+...),... x-°» ( a nn + — + 



x / v x 



where <j\, . . . , u n differ from <r\, . . . , a„ by integers. By a process of 
reduction precisely like that employed earlier (§ 11) one may further 
modify Y (x) and simultaneously Y x (x) so as to preserve all of the 
properties already noted and to finally obtain \a ii \9 i 0. It is to be 
recalled that | Yq (x) | is not identically zero and can at most vanish 
to a finite order at x = oo ; for it is this fact that enables us to conclude 
that the process of reduction terminates. 

The argument of the preceding paragraph shows that Y Q (x) and 
F M (x) will be matrix solutions of a system (54) with coefficients poly- 
nomials in x of degree ju at most. 
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It is deserving of notice that the properties stated above determine 
the location of the poles of the elements of Y x (,r), namely at the zeros 
of | P(x) | within the loop, the maximum order of the pole of any 
element being precisely the order of the corresponding zero of | P(x) \. 
This is an immediate consequence of the relation. 



YJ.r)= Y (x)P- l (z). 
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